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Abstract. - We describe a modified transfer matrix renormalization group (TMRG) 
algorithm and apply it to calculate thermodynamic properties of the one- dimensional 
t-J model. At the supersymmetric point we compare with Bethe ansatz results and 
make direct connection to conformal field theory (CFT). In particular we study the 
crossover from the non-universal high T lattice into the quantum critical regime by 
calculating various correlation lengths and static correlation functions. Finally, the 
existence of a spin-gap phase is confirmed. 



The t-J model is one of the most fundamental systems of strongly correlated 
electrons. The two-dimensional version has attracted much attention because it 
is believed that it describes the basic interactions in the copper-oxygen planes of 
high-Tj. superconductors. For the one-dimensional (ID) t-J model much progress has 
been achieved using various analytical and numerical techniques 
At the supersymmetric point (SUSYP) J /t = 2 the model is solvable by the Bethe 
ansatz (BA) and ground state properties as well as the excitation spectra have been 
obtained exactly ||^ . Because the two critical excitations of spin and charge type are 
separated, they can be described by two independent c—\ Virasoro algebras. By a 
combination of finite-size results from the BA and CFT it is therefore also possible 
to calculate the critical exponents of algebraically decaying correlation functions 
(CF's) confirming Tomonaga-Luttinger liquid (TLL) properties [Q. By numerical 
calculations it was shown that the phase diagram consists of a TLL phase for small 
J /t and a phase separated state for J /t large . In between a spin-gap (Luther- 
Emery) phase was conjectured for low densities |l| and confirmed a few years later 
by different methods ||,^,||. However, the results for the spin-gap phase strongly 
depend on the applied methods and the phase boundaries are controversial. Much 
less is known about thermodynamics of the ID t-J model. Only at the SUSYP 
thermodynamic quantities have been obtained by the BA ||^,^. 

To study thermodynamic properties away from the SUSYP the TMRG provides 
a powerful numerical tool. This method is particularly suited, because the thermo- 
dynamic limit is performed exactly and it has been applied successfully to various 
ID systems before [|o|. The original idea was to decompose a Hamiltonian H with 
nearest neighbour interactions into even (i?e) and odd parts (i?o)- By applying the 
Trotter formula the partition function is expressed as 

Z = Tr c-P" = lim Tr | \c-^"^e'^"^]''] (1) 
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Fig. 1 - The left part shows the usual Trotter mapping of the ID quantum chain to a 
2D classical model with checkerboard structure where the vertical direction corresponds 
to imaginary time. All lattice points of the classical model belong to the physical lattice 
at different imaginary time steps and the QTM is a two-column transfer matrix. The 
right part shows the alternative mapping as described in the text. The classical model 
has alternating rows and additional lattice points in a mathematical auxiliary space. The 
QTM in this formulation is only one column wide. In both figures the shaded plaquettes 
denote the same Boltzmann weight. 



with e — P/M ^ P being the inverse temperature and M an integer Trotter number 
leading to a classical model where the column-to-column transfer matrix (QTM) has 
checkerboard structure (see fig. |^). However, this transfer matrix is unnecessarily 
wide as the repeat length of the classical system is 2. The consequence are several 
disadvantages: (1) The wavcvector A; of a CF is not uniquely determined, i.e. this 
transfer matrix cannot distinguish between k and fc -I- vr; (2) the calculation of 
CF's is complicated, because even and odd as well as distance 1 have to be treated 
separately; (3) the cost of computer memory is unnecessarily large. 

We therefore applied a different Trotter-Suzuki mapping leading to a classical 
lattice with alternating rows (see fig. ^ where the partition function is given by 

lim Trf[ri(6)r2(e)]''''/n with Ti,2(e) = T^^l exp [-ei/ + 0(6^)] (2) 

and Tfi^L being the right- and left-shift operators, respectively. Such a mapping 
is often used in the context of exactly solvable models and is in general applicable 
for ID systems with nearest neighbour interactions. In this formalism the QTM is 
formulated for a single column, thus overcoming the disadvantages of the checker- 
board decomposition. In both transfer matrix approaches the thermodynamic limit 
for a fixed Trotter number M is performed exactly, because the free energy of the 
chain with infinite length is given solely by the largest eigenvalue Aq of the QTM, 
Tm- Within the alternative mapping the free energy is given explicitly by 

/oo,A/ = -rinAo , (3) 

where Ag is unique and a real, positive number for all temperatures. Since a vanish- 
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ing gap between leading and next leading eigenvalue indicates a phase transition, 
such a degeneracy is not possible for a ID quantum system at finite tempera- 
ture. The QTM is a real but non-symmetric matrix with a block structure due to 
conserved quantities (here spin and charge or equivalently the number of up-spin 
electrons and down-spin electrons Ni). The two-point CF is given by 

a 

with Mq, being matrixelements that can be directly evaluated and correlation lengths 
(CL's) and wavevectors ka which are given by nextleading eigenvalues Aq, of the 
QTM. Note that several ^ with the same k appear in the asymptotic expansion. In 
the structure factor each term yields a (measurable) Lorentz function with center 
at ka, height ~ Ma^a/T^ and width ~ The sharpest peak corresponds to the 

leading instability towards the onset of long range order and hence, a crossover in 
the leading CL indicates a change of the nature of the long range order. In addition, 
it is also possible to calculate imaginary time correlations G(r, r) directly within 
the TMRG algorithm ||ll[]. The fatal point is that the analytical continuation is 
an ill-posed problem leading to unreliable results. What is calculated without fun- 
damental problems are static CF's defined by G{r,z = 0) = Jq dr G{r,T). This 
CF can be expressed by the largest eigenvalue and the corresponding left and right 
eigenvectors |*^), l^f^) of the QTM: 

M 

G{r,z = 0) = —^{^^\fMTl-^fM\^^) with fM = Y.TM{A.k) (5) 

for distances r > 1, where TM^A^-k) denotes the usual transfer matrix Tm with 
the considered operator A added at imaginary time position t = e ■ k. The static 
autocorrelation G(r = 0, 2 = 0) has to be treated separately in a similar way. 

The infinite DMRG algorithm is used to increase the Trotter number M being 
equivalent to a decrease of the temperature T. The numerical algorithm is quite 
similar to the one described in Ref. for the checkerboard QTM and will therefore 
not be explained here in detail. In the following, we use a grandcanonical description 
of the t-J model, where the Hamiltonian is given by 

(6) 

with P being the projection operator onto the Hilbert-subspace without double oc- 
cupancy. Note, that the density is a function of chemical potential and temperature 
n = n(/x, T). To check the accuracy of both TMRG algorithms, we have calculated 
several thermodynamic quantities at the SUSYP and compared with BA results. 
With 100 states retained in the DMRG and e = 0.05, the accuracy of the free en- 
ergy in comparison to BA is of the order lO""* if T/t > 0.1 and of the order 10"'^ 
for temperatures down to T/t = 0.01 (see fig. ^) while the deviations between the 
two different algorithms are always an order of magnitude smaller. By applying a 
small magnetic field or a small change in the chemical potential, we have calculated 
spin (xs) and charge susceptibihties (xc) at the SUSYP as shown in fig. ^. Again 
the numerical results are in good agreement with the exact results. Expressing the 
susceptibilities as a sum of two-point CF's, the high temperature asymptotics is 
obtained as Xs ^ 1/6T and Xc ~ 2/9T, also consistent with numerics. CFT yields 
the zero temperature values Xs — l/(27r?;s) and Xc — £,c{QY I [t^Vc) with Vs^c being 
the velocities of the spinon and holon excitations, respectively, and S,c{Q) being 
the dressed charge. We have calculated the velocities as well as the dressed charge 
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Fig. 2 - Deviation of the free energies calculated by means of TMRG in comparison to 
the exact results for different temperatures. The lines are guides to the eye. The inset 
shows the BA results (symbols) plotted on the TMRG results (lines). The three chemical 
potentials correspond to high (/i = —0.7), medium (/i — —1.4) and low — —1.9) electron 
density at low temperatures. 

Fig. 3 - The main figure shows the spin susceptibilities, the inset the charge susceptibilities 
(compressibilities), where again lines denote the numerical results and symbols the exact 
ones. The symbols at T = denote in both graphs the CFT results (see text below). 



by BA leading to the ground-state susceptibilities shown in fig. y. Note, that Vc 
vanishes for the limiting cases n — > or n — > 1 leading to a divergent compress- 
ibility at T = 0. All analytical results together show that the TMRG algorithm 
yields accurate results from the high-temperature limit to temperatures of the order 
T/t = 0.01 corresponding to 2000 renormalization steps. Finally, we want to men- 
tion that neither the described algorithm nor the one based on the checkerboard 
decomposition ever breaks down after performing only a few RG-steps as described 
in Ref. ||l^ although we have varied J/t from the free fermion point to the phase 
separated state and the number of states between 24 and 300. 

The supersymmetric t- J model belongs to the TLL universality class with gapless 
spin and charge excitations meaning there is a quantum critical point at T = 0. 
The linear dispersion of the critical excitations leads to universal low-temperature 
properties like a quadratic temperature dependence of the free energy / ~ bq — aT^, 
a linear regime in the specific heat and CL's diverging as ^ 1/T. More insight has 
been gained from conformal invariance describing the t-J model by two independent 
Virasoro algebras with central charges c — 1 but different Fermi velocities Vc,s- The 
critical exponents of algebraically decaying correlation functions have been obtained 
by the BA due to the relation of finite size energy gaps and scaling operators in 
CFT By the usual conformal mapping of the complex plane onto a cylinder 
these results are easily extended to finite temperatures. However, universality is 
only given in the so called quantum critical regime determined by T <C i. If this 
condition is not fulfilled the properties are non-universal and dependent on the 
microscopic Hamiltonian. To study the crossover from the non-universal high T 
lattice into the quantum critical regime, we focus on the temperature dependence 
of CL's and static CF's. 

Respecting the selection rules, the density-density (d-d) and longitudinal spin- 
spin (s-s) CL's are in the block of the QTM with unchanged quantum numbers 
{ANij — AN^ = 0). To distinguish between them, the matrixelement in eq. (^) 
has to be calculated explicitly. The leading d-d and s-s CL's (times temperature) 
and the corresponding wavevectors for J/t — 2.0 and n = —1.4 {nr^o ~ 0.52) with 
zero magnetic field are shown in fig. ^. The non-universal regime is characterized by 
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Fig. 4 - Temperature dependence of the leading d-d and s-s CL's for J/t — 2.0 and 
H = —1.4. Each curve corresponds to one eigenvalue of the QTM and the respective 
wavevectors may take identical or different values. The triangle up (star) gives the zero 
temperature result from CFT for the non-oscillating s-s (d-d) and the square that for the 
Akp part of the d-d. The triangle down at zero temperature denotes the CFT value for 
the 2kF part of s-s and d-d, whereas the triangles down at T/t = 0.1 are given by CFT 
plus logarithmic corrections as described in the text. The inset shows the wavevectors k 
in the case of incommensurate oscillations. The circles denote the values for k = i:2kF 
at zero temperature as expected from CFT. In the temperature range shown in the inset, 
the density n(p.,T) varies only between 0.51 and 0.52. 

Fig. 5 - Static s-s CF at T = 2.0 (left graph) showing 7r-oscillations and at T = 0.1 
(right graph) showing incommensurate oscillations. The dotted lines denote envelopes 
corresponding to the exponential decay. The straight lines are guides to the eye. 



several erossovers between the CL's driven by temperature. In the high temperature 
hmit all CL's show commensurate oscillations (i.e. k = 0,7r) with a leading tt- 
oscillating s-s CL and a leading non-oscillating d-d CL. In the low-temperature 
regime, T <C t, no crossovers occur and we expect universal properties. Here the 
oscillations of a s-s and a d-d CL are incommensurate with wavevectors depending on 
temperature. Note, that the crossovers to incommensurate oscillations are shifted 
to lower temperatures with increasing particle density and do not occur in the 
Heisenberg limit n = 1. Because the Fermi momentum kF-\{i) for up (down) spin 
electrons at zero temperature is given by fep-m) = 7r(n ± 2m)/2 where m is the 
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magnetization, we identify these CL's as the 2/ci?-oscillating parts. Using the CFT 
results for the critical exponents of the algebraically decaying CF's at T = and 
applying the conformal mapping onto the cylinder, we obtain 



2ttT 



T 



with the scaling dimensions Xc^s- For the non-oscillating part of the s-s (d-d) CL 
we receive 7 — fs/27r (7 — Vc/^tt) whereas 7 = 2-yc/{7r(2^ -I- ^dQY')} for both 
2fcir-oscillating parts. In the d-d correlation there is also a 4A;i;-oscillating part 
with 7 = '!;c/{27r^^((5)}, however, this CL is so small that we have not calculated 
it numerically. For the non-oscillating parts we find a good agreement with the 
numerical results, but the 2/cf s-s and d-d CL's are not equal as expected from CFT. 
This is a consequence of different logarithmic corrections which are not directly 
accessible within BA but known from TLL theory The multiplicative term 

In^^^^ r (In^^^ r) for the 2kp part of the d-d (s-s) correlation at T = can be 
regarded as an effective, distance dependent correction of the scaling dimension x 
at finite temperature 

' x-l'^^^ (a ^-3/2, 1/2), (8) 



2 In 



r 



where the relevant length scale r is the correlation length at the considered temper- 
ature. This correction leads to an excellent agreement with the numerical results 
(see triangles down in fig. ^ . Following only the largest s-s CL a sharp crossover at a 
critical temperature ~ 0.8 from 7r-oscillations to non-oscillations occur. However, 
the matrixelement belonging to the non-oscillating CL is rather small and we expect 
the 2kp CL dominating the asymptotics of the spin CF for T <^Tc. This crossover is 
also visible when the static spin CF is calculated explicitly. At T = 2.0 the CF shows 
TT-oscillations (see left part of fig. ||) and a fit {S^{0)S^{r)) '-^ exp(— r/^) cos(fcr -|- S) 
gives perfect agreement for ^ and k within errors of the order 10^** with the direct 
calculation of CL's. On the other hand, incommensurate oscillations dominate at 
T = 0.1 as shown in the right part of fig. ||. Again the fit values coincide rather 
precisely with the directly calculated values. Although these crossovers at well de- 
fined finite temperatures are visible in physical quantities (i.e. two-point CF's), they 
have nothing to do with phase transitions. We want to point out that any thermo- 
dynamic quantity derived from the free energy is an analytic function at finite T. 
Phase transitions and corresponding singularities only occur at T = 0. However, 
quantities describing the asymptotics of correlation functions (i.e. CL) may show 
non-analyticities even at finite temperature. 

As already mentioned in the introduction a spin-gap or Luther-Emery (LE) 
phase is expected for values of J/ 1 near the phase separated state at least for small 
densities |l|, |[ ^, 0, |^ . Although the TMRG algorithm is not suited to determine 
the precise phase boundaries, because a spin gap A is only visible at temperatures 
T < A, we want to show exemplarily and free of any additional assumptions the 
existence of such a phase. The spin susceptibilities for two different values of J/t 
are shown in fig. ^ where the chosen chemical potentials correspond in both cases 
to n « 0.2 in the low-temperature limit. The quadratic dispersion of a ID gapped 
system leads to x ~ exp{—A/T)/^/T for the low-temperature susceptibihty. Using 
this function for a fit of the numerical data, we find A = 0.05 ± 0.01 in both 
cases. Another proof of LE properties is given by the calculation of s-s and d-d 
CL's (see fig. |7|). In the low-temperature limit the non-oscillating s-s CL seems 
to be finite (T • ^ 0) whereas the non-oscillating d-d CL is not affected and 
diverges as ^ ~ 1/T. However, we are not able to present numerical data for lower 
temperatures, which could support this scenario further. 
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Fig. 7 - Correlation lengths at J/t = 2.9 and n = —2.13. Note, that the oscillations of all 
shown correlation lengths are commensurate over the entire temperature range. 



Summarising, wc have described a modified TMRG algorithm based on a novel 
Trotter-Suzuki mapping where the corresponding QTM is only one column wide. 
The advantages are unambiguously determined wavevectors k, a simplified calcula- 
tion of CF's and a reduction in required computer memory. By calculating thermo- 
dynamic quantities at the SUSYP for different chemical potentials and comparison 
with BA, we have shown that the numerical algorithm yields very accurate results 
down to temperatures of T/t = 0.01 (2000 RG steps). In particular, wc have inves- 
tigated temperature driven crossovers in the correlation lengths and have compared 
the numerical results in the quantum critical regime with CFT where logarithmic 
corrections turned out to be important. For J/t « 3.0 and low densities a spin gap 
was estimated from the spin susceptibilities proving the existence of a LE phase 
what was further supported by the calculation of s-s and d-d CL's. 

* * * 
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